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The authors have recently proposed a "microcanonical functional integral" represen- 
tation of the density of quantum states of the gravitational field. The phase of this 
real-time functional integral is determined by a "microcanonical" or Jacobi action, 
the extrema of which are classical solutions at fixed total energy, not at fixed total 

,— i ' time interval as in Hamilton's action. This approach is fully general but is especially 

O ' 
CO 



well suited to gravitating systems because for them the total energy can be fixed sim- 
ply as a boundary condition on the gravitational field. In this paper we describe how 

CO . 

to obtain Jacobi's action for general relativity. We evaluate it for a certain complex 
metric associated with a rotating black hole and discuss the relation of the result to 
i , the density of states and to the entropy of the black hole. 

1 DEDICATION 

We dedicate this paper to Yvonne Choquet-Bruhat in honor of her retirement fol- 
lowing a brilliant career. JWY would like to thank her for friendship, support, and 
their happy collaboration in "analysis, manifolds, and physics" [1]. 

2 INTRODUCTION 

We are concerned here with the description of a stationary and axisymmetric rotating 
black hole in a closed system in thermodynamic equilibrium with its environment. 
The equilibrating radiation in the closed system surrounding the hole is the "envi- 
ronment". This "radiation fluid" rotates at a constant angular velocity [2] [3]. Unlike 
the more familiar non-relativistic case, here one must take into account explicitly the 
spatial finiteness of the system in order to avoid super-luminal velocities resulting 
from the rotation. Even though we shall ignore the explicit effects of the equilibrat- 
ing radiation in this work, its presence-in-principle must be kept in mind in order 
to give the problem a physically and mathematically reasonable formulation. 



The key conserved quantities in closed systems like ours are the total energy and the 
total angular momentum. Their "Massieu" conjugates [4] are, respectively, inverse 
temperature (3, and (3 multiplied by an appropriate angular velocity. An important 
question is which among these quantities to specify in advance, that is, which "en- 
semble" picture to employ. We choose here to fix energy and angular momentum as in 
the microcanonical ensemble. This choice leads us to employ the general relativistic 
version of Jacobi's action. 

Jacobi's form of the action principle [5] [6] involves variations at fixed energy, rather 
than the variations at fixed time used in Hamilton's principle. The fixed time interval 
in Hamilton's action becomes fixed inverse temperature (3 in the "periodic imaginary 
time" formulation, thus transforming Hamilton's action into the appropriate (imagi- 
nary) phase for a periodic path in computing the canonical partition function from a 
"Euclidean" Feynman functional integral [7] . (We are here and in the next paragraph 
speaking only of energy and inverse temperature in order to simplify the discussion. 
Similar remarks hold for angular momentum and its conjugate.) In contrast, fixed en- 
ergy is suitable for the microcanonical ensemble and, correspondingly, Jacobi's action 
is the phase in an expression for the density of states as a real-time "microcanonical 
functional integral" (MCFI) [8]. 

Let us characterize briefly the canonical and microcanonical pictures. Neither picture 
can hold with perfect precision even in principle when gravity is taken into account 
because the "infinite" heat bath of the canonical picture and the "perfectly adiabatic" 
walls of the microcanonical picture are both at variance with the known physics 
of the gravitational field. However, each picture still provides a useful framework 
for discussion. Furthermore, the long-range, unscreened nature of the gravitational 
interaction means that simple notions of ordinary statistical thermodynamics like 
"extensive", "intensive", etc., no longer apply in general. With gravity, statistical 
mechanics is inherently global. 

In the canonical picture, with a fixed temperature shared by all constituents of a 
system, there are no constraints on the energy. This feature simplifies combinatorial 
(counting) problems and leads to factorization of the partition function for weakly 
coupled constituents. For gravitating systems in equilibrium, the temperature is 
not spatially uniform because of red-shift and blue-shift effects. In such cases, the 
relevant temperature is that determined at the spatial two-boundary B of the system 
[9]. It can be specified by a boundary condition on the metric [10] [11]. It is then 
used in conjunction with Hamilton's principle, which is the form of the gravity action 
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in which the metric is fixed on the history of the spatial two-boundary [12]. (The 
metric determines the lapse of proper time along the history of the boundary.) The 
case of a rotating, charged, stationary, axisymmetric black hole has been treated with 
(grand) canonical boundary conditions in [13]. See also [10] and [11]. 

With its constraint on the energy, the microcanonical picture leads to stability prop- 
erties more robust than in the canonical case. However, the energy constraint can 
complicate combinatorial problems because the constituents of the system must all 
share from a common fixed pool of energy. For field theories, with a continuous in- 
finity of degrees of freedom, the energy constraint restricts the entire phase space of 
the system unless gravity is taken into account. For gravitating systems, as a conse- 
quence of the equivalence principle, the total energy in a given finite region, including 
that of matter fields, can be given as an integral of certain geometrically well-defined 
derivatives of the metric over a finite two-surface bounding the system. In other 
words, one can find and employ a suitable expression for "quasi-local" gravitational 
energy that does not require an appeal to asymptotic regions. This quasi-local en- 
ergy has the value, per unit proper time, of the Hamiltonian of the spatially bounded 
region, as discussed in detail in [12]. Therefore, if we specify as a boundary condition 
the energy per unit two-surface area, we have constrained the total energy simply by 
a boundary condition [8] [12]. Thus, through the mediation of the gravitational field, 
the canonical and microcanonical cases are placed on similar footing and differ only 
in which of the conjugate variables [14], inverse temperature or energy, is specified 
on the boundary. The corresponding functional integrals, for the partition function 
or density of states, differ in which action gives the correct phase, Hamilton's or 
Jacobi's. We regard the MCFI for the density of states as the more fundamental [8]. 

We shall now outline a recent application of the above reasoning to the case of an 
axisymmetric stationary black hole. The MCFI, in a steepest descents approxima- 
tion, shows that the density of states is the exponential of one-fourth of the area 
of the event horizon, thus confirming in this approximation the Bekenstein-Hawking 
expression for the entropy of a black hole [8]. Full details of the following are given 
in [8] and [12]. (An application of the MCFI to a simple harmonic oscillator-without 
gravity-has been given recently by the authors [15]. No approximations were required 
in this case and the exact energy spectrum was obtained.) 

3 JACOBI'S ACTION FOR GRAVITY 

We here analyze the action for pure general relativity. The method, however, can be 
just as well applied when matter is included and/or for higher-derivative theories of 
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gravity. 



Consider a region of spacetime M = E x I where E are spatial slices and I is an 
interval of the real line. The two-boundary of space E is denoted by B, whose history 
is 3 B = B x /. The orthogonal intersection of any E with 3 B is the two-boundary 
B at some time. Thus a generic B can be regarded as being embedded in E with a 
spacelike unit normal tangent to E at B, and also as being embedded in 3 B with 
a timelike unit normal such that u^n^ = 0. The subspaces of M that correspond 
to the endpoints of / are spacelike hypersurfaces t = t' and t = t" . The notation 
f t , d 3 x denotes an integral over t" minus an integral over t'. The spacetime metric 
is and its scalar curvature is 3?. Then "Hamilton's" action for general relativity, 
in which the metric is fixed on the boundary, is given by a variant of the Hilbert 
action, namely [12] [16] 

t" 

S[g] = J- I d 4 x^(^t - 2A) + - / d 3 xVhK - - I d 3 x^j® - S° , (1) 
2k Jm k Jf K Jm 

where the metric and extrinsic curvature of the slices E are hij and while those 
for 3 B are 7^ and 0^ . (Latin letters i, j, ... are used as tensor indices for both E 
and 3 B. No cause for confusion arises from this convention.) In (1), k = 87tG and we 
set Newton's constant G = 1 henceforth. The term S° in (1) is a functional of the 
metric 7^ of 3 B. The purpose of this term is to determine the "zero" of energy and 
momentum. It will not turn out to affect the Jacobi action for general relativity. (S° 
is analyzed in [12] but was not included in [16].) 

The canonical form of "Hamilton's" action (1) is [12] 

S=[ d 4 x[P lJ h tJ - NH-V l Hi] - [ d 3 x^[Ne - V 1 ^] , (2) 

where P* J is the Arnowitt-Deser-Misner momentum conjugate to h^, and Ti and Tii 
are the usual Hamiltonian and momentum constraints. The lapse function is denoted 
by and the shift vector by V 1 . In the surface term of (2), a denotes the determinant 
of the two-metric cr^ induced on B as a surface embedded in E; likewise, ni denotes 
B's unit normal in E and kij denotes the corresponding extrinsic curvature. The 
energy surface density e and momentum surface density ji are given by [12] 



The total quasi-local energy is the integral of e over B. In obtaining (2), we have 
assumed that S°, if present, is a linear functional of N and V 1 on 3 B. In the total 
Hamiltonian extracted directly from (2), the shift vector must satisfy riiV' l \B = 0. 
Each of these points is discussed in [12]. 

In the action (2), one varies hij, , N, and V 1 to obtain 

t" 

5S = (terms giving the equations of motion) + / d 3 x P^Shij 

Jt> 

-[ d 3 x v ^\e5N-jJV a -(N/2)s ab 5a ab \ , (4) 

where indices a, 6, ... are used to denote tensors on B, i.e., E-tensors that are 
orthogonal to n l . The surface stress tensor s ab on B is given in [12] and does not 
concern us here. We see from (4) that suitable boundary conditions for S are obtained 
by fixing the metric induced on the boundary elements t', t" , and S B of M. In 
particular, the lapse function N is fixed on 3 B, where it determines proper time 
elements Ndt on 3 B along the unit normal associated with the foliation of 3 B by 
B; hence, we refer to S as "Hamilton's action" in canonical form. 

What we define as the microcanonical action S m is, in essence, Jacobi's action for 
general relativity. It is obtained from S by a canonical transformation that changes 
the appropriate boundary conditions on 3 B from fixed metric components N, V a , and 
a a h to fixed energy surface density s, momentum surface density j a , and boundary 
metric o ao . Thus, define [8] 



(5) 



S m = S+ [ d 3 x v ^[Ne-V a j a ] 

J 3 B 

= [ d^x [P lJ h l3 — NH — V l Hi] . 
Jm 

From (4) , it follows that the variation of S m is 

t" 

SS m = (terms giving the equations of motion) + / d 3 xP^dhij 

Jv 

+ [ d 3 x\N5(V^e)-V a 5(V^ja) + (NV^/2)s ab 5a ab \ . (6) 

J 3 B 1 J 

This result shows that solutions of the equations of motion extremize S m with respect 
to variations in which e, j a , and a ao are held fixed on the boundary B. Observe 
that the unspecified subtraction term S° does not appear in S m . Nevertheless, the 
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variation (6) of S m involves e, j a , and s ab , which do depend on S° through their 
definitions. However, all dependences on S° in the boundary variation terms of SS m 
actually cancel because of the requirement that 5° be a linear functional of the lapse 
and shift [8] [12]. Thus, neither S m nor its variation depends on S°. In other words, 
as long as energy and momentum are to be fixed, as in S m , a quantity like S° whose 
only role is to determine their "zero points" is irrelevant. 

4 MICROCANONICAL FUNCTIONAL INTEGRAL 

In [8] and [15] we showed that for nonrelativistic mechanics the density of states 
is given by a sum over periodic, real time histories, where each history contributes 
a phase determined by Jacobi's action. In the case of nonrelativistic mechanics, 
the energy that is fixed in Jacobi's action is just the value of the Hamiltonian that 
generates unit time translations. For a self-gravitating system, the Hamiltonian has a 
"many-fingered" character: space can be pushed into the future in a variety of ways, 
governed by different choices of lapse function N and shift vector V 1 . The value 
of the Hamiltonian incorporated into (2) depends on this choice. More precisely, 
the value of the Hamiltonian is determined by the choice of lapse and shift on the 
boundary B, since the lapse and shift on the domain of E interior to B are Lagrange 
multipliers for the (vanishing) Hamiltonian and momentum constraints. Accordingly, 
the energy surface-density e and momentum surface-density j a for a self-gravitating 
system play a role that is analogous to energy for a nonrelativistic mechanical system. 

The above considerations lead us to propose that the density of states for a spatially 
finite, self-gravitating system is a functional of the energy surface-density e and mo- 
mentum surface-density j a . In addition to these energy-like quantities, the density 
of states is also a functional of the metric a a b on the boundary B, which specifies 
the size and shape of the system. In the absence of matter fields, these make up the 
complete set of variables and v[e, j a , cr &] is interpreted as the density of quantum 
states of the gravitational field with energy density, momentum density, and bound- 
ary metric having the values e, j a , and a a b- The action to be used in the functional 
integral representation of v is S m , which describes the gravitational field with fixed 
e, j a , and a a b. Note that e, j a , and a a b replace the traditional thermo dynamical ex- 
tensive variables. Our variables are all constructed from the dynamical phase space 
variables (hij, P iJ ) for the system, where the phase space structure is defined using 
the foliation of M into spacelike hypersurfaces. (We expect this to be a defining 
feature of extensive variables for general systems of gravitational and matter fields.) 

We propose [8] that the density of states of the gravitational field is defined formally 
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by 

u[e,j,a] =J2[ VHexp{iSm) . (7) 

M ^ 

(Planck's constant has been set to unity.) The sum over M refers to a sum over 
manifolds of different topologies. The three-boundary for each M is required to have 
topology dM = B x S . If B has two-sphere topology, then the sum over topologies 
includes M = (ball) x S 1 , with dM = d(ball) x S 1 = S 2 x S 1 . Another example is 
M = (disk) x S 2 , with dM = <9(disk) x S 2 = S 1 x S 2 . The action S m that appears 
in Eq. (7) is the microcanonical action (5) of the previous section applied to the 
manifolds M with a single boundary component dM = 3 B = B x S 1 . The functional 
integral (7) for v is a sum over Lorentzian metrics g^ v . Note that the microcanonical 
action may require the addition of terms that depend on the topology of M, such as 
the Euler number. 

In the boundary conditions on dM = B x S 1 , the two-metric a a b that is fixed on 
the hypersurfaces B is typically real and spacelike. Likewise, the energy density 
e is real, which requires the unit normal to dM to be spacelike. Therefore, the 
Lorentzian metrics on M must induce a Lorentzian metric on dM, where the timelike 
direction coincides with the periodically identified S . Note, however, that there are 
no nondegenerate Lorentzian metrics on a manifold with topology M = (disk) x S 2 
that also induce such a Lorentzian metric on dM. This implies that the formal 
functional integral (7) for the density of states must include degenerate metrics. (For 
a discussion of the role of degenerate metrics in classical and quantum gravity, see 

[17]-) 

Now consider the evaluation of the functional integral (7) for fixed boundary data 
£ 5 jai °~ab that correspond to a stationary, axisymmetric black hole. That is, start 
with a real Lorentzian, stationary, axisymmetric, black hole solution of the Einstein 
equations, and let T = constant be stationary time slices that contain the closed 
orbits of the axial Killing vector field. Next, choose a topologically spherical two- 
surface B that contains the orbits of the axial Killing vector field, and is contained 
in a T = constant hypersurface. From this surface B embedded in a T = constant 
slice, obtain the data e, j a , and a a b- In the functional integral for u[e,j,a], fix this 
data on each t = constant slice of dM. Observe that, to the extent that the physical 
system can be approximated by a single classical configuration, that configuration 
will be the real stationary black hole that is used to induce the boundary data. 

The functional integral (7) can be evaluated semiclassically by searching for four- 
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metrics that extremize S m and satisfy the specified boundary conditions. Ob- 
serve that the Lorentzian black hole geometry that was used to motivate the choice of 
boundary conditions is not an extremum of S m , because it has the topology [Wheeler 
(spatial) wormhole] x [time] and cannot be placed on a manifold M with a single 
boundary S 2 x S 1 . However, there is a related complex four-metric that does ex- 
tremize Sm, and is described as follows. Let the Lorentzian black hole be given 
by 

ds 2 = —N 2 dT 2 + h l3 \dx l + V l dT)(dx J + V J dT) , (8) 

where N, V 1 , and hij are T-independent functions of the spatial coordinates x l . The 
horizon coincides with N = 0. Choose spatial coordinates that are "co-rotating" with 
the horizon. Then the proper spatial velocity of the spatial coordinate system relative 
to observers at rest in the T = constant slices vanishes on the horizon, (V l /N) = 0, 
and the Killing vector field d/dT coincides with the null generator of the horizon. 
As shown in [13], the complex metric 

ds 2 = -(-iNfdT 2 + h l3 (dx l - iV l dT)(dx j - iV J dT) , (9) 

where the coordinate T is real, satisfies the Einstein equations everywhere on a 
manifold with topology M = (disk) x S 2 , with the possible exception of the points 
N = where the foliation T = constant degenerates. The locus of those points 
N = is a two-surface called the "bolt" [18]. Near the bolt, the metric becomes 

ds 2 w N 2 dT 2 + h tj dx i dx j , (10) 

and describes a Euclidean geometry. The sourceless Einstein equations are not sat- 
isfied at the bolt if this geometry has a conical singularity in the two-dimensional 
submanifold that contains the unit normals h l to the bolt for each of the T = constant 
hyper surfaces. However, there is no conical singularity if the circumferences of circles 
surrounding the bolt initially increase as 2ir times proper radius. The circumference 
of such circles is given by PN, where P is the period in coordinate time T. There- 
fore the absence of conical singularities is insured if the condition P(h l diN) = 2n 
holds at each point on the bolt, where h 1 is the unit normal to the bolt in one of 
the T = constant surfaces. Because the unit normal is proportional to diN at the 
bolt, this condition restricts the period in coordinate time T to be P = 2ty/k h , where 
k h = [{diN)h 13 (djN)] 1 ^ 2 \ H is the surface gravity at the horizon of the Lorentzian 
black hole (8). 

The lapse function and shift vector for the complex metric (9) are iV = — iN and 
V 1 = —iV 1 . Thus, (9) and the Lorentzian metric (8) differ only by a factor of — i 
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in their lapse functions and shift vectors. In particular, the three-metric hij and its 
conjugate momentum P lj coincide for the stationary metrics (8) and (9) [13]. Since 
the boundary data s, j a , and a a b are constructed from the canonical variables only, 
the complex metric (9) satisfies the boundary conditions imposed on the functional 
integral for v[e, j, a] . 

The complex metric (9) with the periodic identification given above extremizes the 
action S m and satisfies the chosen boundary conditions for the density of states 
v[s,j, a]. Although this metric is not included in the sum over Lorentzian geometries 
(7), it can be used for a steepest descents approximation to the functional integral 
by distorting the contours of integration for the lapse N and shift V % in the complex 
plane. Then the density of states is approximated by 

v[e,j,a] wesxp(iS m [-iN,-iV,h]) , (11) 

where S m [— iN, — iV, h] is the microcanonical action evaluated at the complex ex- 
tremum (9). The density of states can be expressed approximately as 

v[e,j,a] « exp(S[e,j,o]) , (12) 

where S[e,j,a] is the entropy of the system. Then the result (11) shows that the 
entropy is 

S[e,j,a}^iS m [-iN,-iV,h] (13) 
for the gravitational field with microcanonical boundary conditions. 

In order to evaluate S m for the metric (9), we start with the microcanonical action 
written in spacetime covariant form [8] and perform a canonical decomposition under 
the assumption that the manifold M has the topology of a punctured disk xS 2 . That 
is, the spacelike hypersurfaces £ have topology I x S 2 , and the timelike direction is 
periodically identified (S 1 ). The outer boundary of the disk corresponds to the three- 
boundary element 3 B of M (denoted dM previously) on which the boundary values 
of e, j a , and a ao are imposed. The inner boundary of the disk, the boundary of the 
puncture, appears as another boundary element S H for M. (No data are specified at 
3 H.) The canonical decomposition results in [8] 

S m = I d 4 x[P lJ h lJ -NH-V l H i \+ [ d 3 xV^[n l (d l N)/ K + 2n t V J P lJ /Vh] , (14) 
Jm Jm 

where the expression ai = (diN)/N for the acceleration of the timelike unit normal 
has been used. The boundary term at 3 H was first given in [13]. 
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Now evaluate the action S m on the punctured disk xS 2 for the complex metric 
(9), and take the limit as the puncture disappears to obtain a manifold topology 
M = (disk) x S 2 . In this limit, the smoothness of the complex geometry is assured 
by the regularity condition on the period of T. Since the metric satisfies the Einstein 
equations, the Hamiltonian and momentum constraints vanish, and the terms P l3 hij 
also vanish by stationarity. Moreover, the second boundary term at 3 H is zero because 
the shift vector vanishes at the horizon. Thus, only the first of the boundary terms 
at 3 H survives. Evaluating this term for the complex metric (9), that is, for the lapse 
function N = — iN, and using the regularity condition P = 2tv/k h , we find for the 
micro canonical action 

S m [-iN, -iV, h] = -- f dT [ d 2 xV&n l d z N = --A H , (15) 

k Jo J 4 

where A H is the area of the event horizon for the Lorentzian black hole (8). 

The result (15) for the microcanonical action evaluated at the extremum (9) leads to 
the following approximation for the entropy (13): 

S[e,j,a]*^A H . (16) 

The generality of the result (16) should be emphasized: The boundary data s, j a , 
and a a b were chosen from a general stationary, axisymmetric black hole that solves 
the vacuum Einstein equations within a spatial region with boundary B. Outside 
the boundary £?, the black hole spacetime need not be free of matter or be asymp- 
totically flat. Thus, for example, the black hole can be distorted relative to the 
standard Kerr family. Furthermore, recall that the quantum-statistical system with 
this boundary data is classically approximated by the physical black hole solution 
that matches that boundary data. The result (16) shows that the entropy of the 
system is approximately 1/4 the area of the event horizon of the physical black hole 
configuration that classically approximates the contents of the system. It also should 
be emphasized that the entropy is independent of the term S° in (1) as has been 
shown in the framework of the canonical partition function [9] . 

Expresson (16) is the principal result we wished to demonstrate here. The physical 
and mathematical limitations of our analysis and possible ways to overcome them are 
described in [8] , which also discusses canonical and grand canonical boundary condi- 
tions for the rotating black hole. A further elaboration of the physical underpinnings 
and implications of our analysis of relativistic rotating systems is given in [19]. 



10 



The authors gratefully acknowledge research support received from the National Sci- 
ence Foundation, grant number PHY-8908741. 

REFERENCES 

[1] Y. Choquet-Bruhat and J.W. York, The Cauchy Problem, pp. 99-172 in General 
Relativity and Gravitation, edited by A. Held (Plenum Press, New York, 1980). 

[2] W. Israel and J.M. Stewart, in General Relativity and Gravitation. II, edited 
by A. Held (Plenum Press, New York, 1980). 

[3] V. Frolov and K.S. Thorne, Phys. Rev. D39, 2125 (1989). 

[4] H.B. Callen, Thermodynamics (Wiley, New York, 1985). 

[5] C. Lanczos, The Variational Principles of Mechanics (University of Toronto 
Press, Toronto, 1970). 

[6] J.D. Brown and J.W. York, Phys. Rev. D40, 3312 (1989). 

[7] R.P. Feynman and A.R. Hibbs, Quantum Mechanics and Path Integrals (Mc- 
Graw-Hill, New York, 1965). 

[8] J.D. Brown and J.W. York, Phys. Rev. D47, 1420 (1993). 

[9] J.W. York, Phys. Rev. D33, 2092 (1986). 

[10] B.F. Whiting and J.W. York, Phys. Rev. Lett. 61, 1336 (1988). 

[11] H.W. Braden, J.D. Brown, B.F. Whiting, and J.W. York, Phys. Rev. D42, 
3376 (1990). 

[12] J.D. Brown and J.W. York, Phys. Rev. D47, 1407 (1993). 

[13] J.D. Brown, E.A. Martinez, and J.W. York, Phys. Rev. Lett. 66, 2281 (1991). 

[14] J.D. Brown, G.L. Comer, E.A. Martinez, J. Melmed, B.F. Whiting, and J.W. 
York, Class. Quantum Grav. 7, 1433 (1990). 

[15] J.D. Brown and J.W. York, "Jacobi's action and the density of states" , to appear 

11 



in Festschrift for Dieter Brill, edited by B.L. Hu and T. Jacobson (Cambridge 
University Press, Cambridge, 1993). 

[16] J.W. York, Found. Phys. 16, 249 (1986). 

[17] G.T. Horowitz, Class. Quantum Grav. 8, 587 (1991). 

[18] G.W. Gibbons and S.W. Hawking, Commun. Math. Phys. 66, 291 (1979). 

[19] J.D. Brown and J.W. York, "Hamiltonian and boundary conditions for rotating 
general-relativistic equilibrium systems" , to be published. 



12 



